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ABSTRACT 

We conduct a detailed leasing, dynamics and stellar population analysis of nine mas¬ 
sive lens early-type galaxies (ETGs) from the X-Shooter Lens Survey (XLENS). Com¬ 
bining gravitational lensing constraints from HST imaging with spatially-resolved 
kinematics and line-indices constraints from VLT X-Shooter (XSH) spectra, we in¬ 
fer the low-mass slope and the low cut-off mass of the stellar Initial Mass Eunction 
(IME): X250 = 2 . 37 lo!i 2 Mow, 250 = 0 - 1311 o!o 26 ^©5 respectively, for a reference 

point with = 250kms“^ and i?eff= 10 kpc. All the XLENS systems are consistent 
with an IME slope steeper than Milky Way-like. We find no significant correlations 
between IME slope and any other quantity, except for an anti-correlation between 
total dynamical mass density and low-mass IME slope at the 87% CL [dx/d\og{p) = 
—O.lOtoiis]* This anti-correlation is consistent with the low redshift lenses found by 
Smith et al. (2015) that have high velocity dispersions and high stellar mass densities 
but surprisingly shallow IME slopes. 

Key words: dark matter - galaxies: ellipticals and lenticular, cD - gravitational 
lensing: strong - galaxies: kinematics and dynamics - galaxies: structure - galaxies: 
formation 


1 INTRODUCTION 


Low mass stars (M ^ 0.3 M©) contribute at most 5-10 per 
cent to the total optical light of the integrated spectrum 
of early-type galaxies (ETGs) but dominate the total stel¬ 
lar mass budget ( Wort he;^| 1994 ). Variations in the number 
of M-dwarfs can therefore dramatically change the stellar 
mass-to-light ratio (T*) without signihcantly affecting the 
characteristics of the observed light (e.g., Bell & de Jong 


2001 Conroy van Dokkum| 2012a Conroy 2013). Hence 
the initial mass function (IMF), and especially its low-mass 
slope, plays a fundamental role in determining an ETC’s 
T* and is crucial to understand the internal structure, the 
formation times and the evolution of these massive galaxies 


(e.g. Blumenthal et al. 1984), providing important insight 


into galaxy evolution mechanisms (e.g. Padmanabhan et al. 


|2004| [Tortora et al.|2009| [Auger et al.|2010| [Graves &: Faber 


* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile. 
(P086.A-0312, PI: Koopmans and P089.A-0364, PI: Spiniello) 
t E-mail:spini@mpa-garching.mpg.de 


2010 Barnabe et al.|20TT ). 


Recent observations, based on gravitational lensing, dy¬ 
namics and/or simple stellar population (SSP) modeling of 
galaxy spectra, indicate that the number of low-mass stars 
in the central region of ETGs increases rapidly with stellar 


velocity dispersion (Treu et al.||2010 

van Dokkum V Con- 

roy||2010| ISpiniello et al.||2011| |2012| 

Cappehari et al.||20121 

La Barbera et al.||2013| 

Tortora, Romanowsky V Napoli-] 

tano|2013 

Spiniello et al. 

2014 Martm-Navarro et al.|2015). 


Recently, however, Smith, Lucey V Conroy (2015, hereafter 
S15) presented results on three very massive lens ETGs with 
high velocity dispersions (cr > 300 kms“^) and a-enhanced 
abundances whose T* are compatible with a Milky Way-like 
IMF slope (although a Salpeter IMF cannot be excluded). 
The potential tension between different studies could indi¬ 
cate that stellar mass (or velocity dispersion) is not the main 
driver behind IMF variations, but only a proxy for other, 
more physically motivated causes. 

Strong gravitational lensing combined with stellar dy¬ 
namics (GL+D) allows one to measure stellar masses of 
galaxies independently from assumptions about their stel¬ 
lar IMF or knowledge of their stellar population. Therefore, 
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Table 1. Properties of the nine XLENS systems used in the analysis. 


XLENS System 

Zlens 

Zsource 

Reff(kpc) 

REin (kpc) 

Lv {W^^LqY 

crxsH (km/s)^ 

IMF slope^ 

SDSSJ0037-0942 

0.1955 

0.6322 

7.03 

4.95 

1.09 + 0.06 

277 + 6 

2.5 + 0.25 

SDSSJ0044+0113 

0.1196 

0.1965 

5.56 

1.72 

0.61 + 0.05 

260 + 8 

- 

SDSSJ0216-0813 

0.3317 

0.5235 

12.6 

5.53 

1.82 + 0.05 

327 + 19 

2.9 + 0.35 

SDSSJ0912+0029 

0.1642 

0.3239 

10.8 

4.58 

1.47 + 0.05 

325 + 10 

2.6 + 0.30 

SDSSJ0935-0003 

0.3475 

0.4670 

20.7 

4.26 

2.07 + 0.07 

380 + 22 

- 

SDSSJ0936+0913 

0.1897 

0.5880 

6.61 

3.45 

0.83 + 0.05 

256 + 18 

— 

SDSSJ0946+1006 

0.2219 

0.6085 

8.33 

4.95 

0.66 + 0.06 

300 + 22 

2.1 + 0.15 

SDSSJ1143-0144 

0.1060 

0.4019 

9.21 

3.27 

0.95 + 0.05 

287 + 18 

2.4 + 0.18 

SDSSJ1627-0053 

0.2076 

0.5241 

6.66 

4.18 

0.79 + 0.05 

303 + 23 

2.3 + 0.26 


1 ) Based on WFPC2 F606W observations. 2 ) Luminosity-weighted stellar kinematics of the lens extracted from a rectangular 
aperture of 2 " x 1^5 centered on the galaxy’s center. See Paper II. 3 ) Best-fit IMF slope from spectroscopic, line-index based 
SSP analysis. For three systems, the spectra do not have S/N high enough to securely constrain the IMF. 


when combined with a measurement of the galaxy luminos¬ 
ity, one obtains a measure of the value of T* which is a func¬ 
tion of the IMF, allowing the latter to be constrained (e.g. 


Treu et ar]|2010| |Spiniello et al.pOllI [Barnabe et al.|[2012| 


2013). Applications of GL+D have shown that the total T 
of massive ETGs increases monotonically with the velocity 
dispersion of the galax y (e.g.|Grillo et al.|2009 Auger et al. 


2010||Treu et al.|2010||Barnabe et al.|2011| [Dutton, Mendel 


& Simard f2012 ). A steepening of the low-mass end of the 

stellar IMF with galaxy mass could therefore be largely re¬ 
sponsible for this change in T*, although a change in the 
dark matter (DM) fraction could also play a role (Auger 
et al.||2009| . 

With the X-Shooter Lens Survey (XLENS: |Spiniello 
et aL|[20TT Barnabe et al.|[2013 hereafter B13) we aim to 
probe variations of the low-mass end of the IME and disen¬ 
tangle stellar and DM in the internal region of lens ETGs, 
via a joint lensing+dynamics+stellar population analysis. 

Except for B13, all studies to date have assumed that 
the lower mass limit of the IME is constant, since varying 
this parameter does not strongly impact the spectra nor 
change the line-index measurements for any assumed IME 
slope. Stars with masses below ~ 0.15 M© are however crit¬ 
ical to determine T*. The only possible way to break this 
degeneracy is by obtaining a measurement of the total stel¬ 
lar mass from an independent method, which, as in B13, is 
provided by a combined lensing and dynamics analysis. 

In this Letter we study the low-mass end properties of 
the stellar IME, inferring its slope and its low-mass cut-off 
(Miow) in nine massive lens ETGs, quadrupling the sample 
presented in B13, which we subsequently augment by three 
recently discovered low-redshift lenses (S15). These two sets 
allow us, for the first time, to study trends of the IME pa¬ 
rameters (slope and low-mass limit) with stellar velocity dis¬ 
persion, effective radius and galaxy mass density. 


2 DATA SET AND ANALYSIS 

The XLENS sample consists of eleven massive lens ETGs 
with cr* = 250 — 400 kms“^, all with multi-band HST and 
VLT XSH data. Here we use nine out of eleven systems 
with the necessary total luminosity and high signal-to-noise 
(S/N) spectra. Gharacteristics of the systems used in this 
work are presented in Table 1. 


Eor the GL+D analysis we make use of the high- 
resolution F 8 I 4 W images of the lens systems, originally se¬ 


lected from the Sloan Lens AGS Survey (SLAGS, Bolton 
et al.|[2008 ). Eor the SSP modeling, UVB-VIS high signal- 


to-noise spectra (S/N > 50) have been obtained during two 
runs between 2011 and 2013 in slit mode (UVB: R = 3300 
with F.'6 X 11'.' slit; VIS: R = 5400, with l'.'5 x 11'.' slit). 

We refer to Paper II (Spiniello et al. 2015, submitted) 
for a more detailed presentation of the full survey sample and 
the spatially resolved kinematics analysis of the systems, to 
B13 for a detailed description of the precise methodology 
and to Paper III (Barnabe et ah, in prep.) for a more ex¬ 
tensive discussion and results of the GL+D analysis of the 
full sample as well as detailed constraints on the DM halo 
fraction and density profiles. 

To infer the IME slope and cut-off mass, we apply the 
same methodology introduced in B13, summarized below. 

i) Erom a GL+D analysis, we determine the posterior 
probability density function (PDE) of the stellar mass of 
each of the nine lens galaxies. This measurement is inde¬ 
pendent from any SSP or photometric analysis and does 
not make any assumption about the IME. As in B13, we use 


the fully Bayesian GAULDRON code (detailed in Barnabe 


& Koopmans 


2007 


and Barnabe et al.| 2012) which is de¬ 
signed to conduct a self-consistent combined modelling of 
both the lensing and kinematic constraints. In particular we 
use a flexible two-component axially-symmetric mass model 
for each lens. We use a generalized Navarro-Erenk-White 
(gNFW) for the DM halo profil(0 consistent with cosmo¬ 
logical simulation findings, as well as observations, and we 
infer the density profile of the luminous mass component 
by deprojecting the multi-Gaussian expansion (MGE, |Em- 


sellem, Monnet & Bacon 1994 Gappellari 2002) fit to the 


observed surface brightness distribution of the galaxy. The 
total stellar mass then sets the normalization of the stellar 
mass distribution and is used to simultaneously model both 
the lensing data set and the stellar kinematic observables. 

ii) We infer a value of the low-mass IME slope from 
XSH spectra, using an extended version of the SSP models 


^ The inner slope of the dark halo density profile is allowed to 
vary between 0 and 2. This accounts for a wide variety of possible 
DM inner slopes, from the pure NFW profile, to both steeper 
(‘contracted’) and shallower (‘cored’) profiles, proposed to try to 
include the effect of baryons on the dark halo inner regions. 
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Figure 1. Marginalized posteriors for the model fitting of IMF slope (x) and cut-off mass as function of a and J?eff (ia unit of 

250 kms“^ and 10 kpc, respectively), including the XLENS systems only (left, red contours) or the XLENS+S15 systems (right, blue 
contours). 


of Conroy & van Dokkum (2012a), as described in Spiniello 


et al. ( 2014| and* Spiniello, Trager & Koopmans (2015). This 

inference (see Table 1) is used as Gaussian prior in our sub¬ 
sequent analysis. Only six out of nine ETGs have spectra 
that provide a solid measurement. For the three remaining 
systems, we assume a flat prior on the IMF slope between 
X = 1.8 and 3.5. 

iii) We use a grid of T* values from the Dartmouth 


Stellar Evolution Program (Chaboyer et al. 2001) as func¬ 
tions of the IMF slope {x = l.S — 3.5) and low-mass limit 
(Miow= 0.01 — 0.22). We use a Gaussian probability dis¬ 
tribution function for total V-band luminosities (Ly) of 
the lens galaxies, as determined from HST data by |Auger| 
et al. (2009), to convert mass-to-light ratios (T*) into stellar 


masses as function of the IMF slope and low-mass cut-off. 

iv) We subsequently use a Markov Chain Monte Carlo 
(MCMC) analysis of the model parameters to derive the 
GL+D+SSP posterior PDFs of the IMF slope and low-mass 
cut-off for each lens galaxy. In the model, x, Miow and Ly all 
vary. The GL+D PDFs for the stellar mass produce the like¬ 
lihood, which we then multiply with the priors on luminosity 
and IMF slope from the HST and XSH data, respectively. 
Finally, the luminosity of the galaxy is marginalized out. 

In this way, we obtain a posterior probability density 
for each system: 

P(x, Miow|data) = J £(dataGL+D|M* = Ty{x,Miow) 

X Lv) • Phst(Lv) • Pssp{x) • P{Miow)dLy (1) 


Tfiow — Tfiow,250 T ^iow,i logj^Q(cr25o) T 2 

(3) 

where X 250 and Miow ,250 represent the slope and the cut¬ 
off mass of the IMF respectively, and (7250 and Rio are the 
velocity dispersion in unit of 250 kms“^ and the effective 
radius in units of 10 kpc. 

We then sample, using a MCMC Metropolis algorithm, 
ah six free parameters (the two biases and four trends with 
velocity dispersion and effective radius) and create their re¬ 
spective univariate and bivariate marginalised posteriors, as¬ 
suming hat priors. The likelihood value is the product of 
the nine LHs values 11^=1 9 P{^^ Miowldata^) times the pri¬ 

ors on the model parameters for the trend function inferred 
above (ah assumed to be hat). The small measurements er¬ 
rors on the effective radius, velocity dispersion and redshift 
are negligible in this process. Each posterior is built from 
10^ samples. 

Besides the nine XLENS galaxies, we augmented the 
MCMC analysis with three high-velocity dispersion, low- 
redshift lenses presented in SI5. Using their constraints on 
the T* and the ‘mass excess factor’ a (see Table 3 in SI5), we 
infer the corresponding value of the IMF slope (unimodal, 
single power-law) for each galaxy. This produces a tightening 
of most of the marginalised posteriors, but there is no major 
shift in the results, indicating that the model is consistent 
with both data-sets and their combination. 


3 MCMC RESULTS 


With these nine P(x, Miowldata^) in hand, which are kept 
hxed in the subsequent analysis, we proceed with inferring 
whether there are any trends of the IMF parameters with 
stellar velocity dispersion (cr) and effective radius (Peff) or 
the total (dynamical) density p = We hrst construct 

a multivariate model: 

X = X250 + x’l logio(o-25o) + X^ logio(Plo) (2) 


Figure shows the univariate and bivariate marginalized 
posteriors resulting from the MCMC model htting of the 
nine XLENS (left, red plots) and of the XLENS+S15 sys¬ 
tems (right, blue plots). The upper panels of the hrst two 
columns show the bias values of the slope and the low- 
mass cut-off dehned for a reference system with a = 250 
kms“^ and Peff = 10 kpc. We hnd X 250 = 2.411q‘J 5 and 
Tflow,250 = 0.1441° “4 Mo for XLENS and *250 = 
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Figure 2. 2-D marginalized PDFs for the model fitting of IMF slope (x) and low-mass cut-off as function of total density 

(p = for the XLENS galaxies only (left, red contours) and for the XLENS+S15 sample (right, blue contours). The data suggest 

an anti-correlation between IMF slope and density. 


and Miow. 250 = O.13llo;o26A^0 for the XLENS+S15 case, 
respectively. The results are fully consistent. In both cases 
a Milky Way-like IMF (x = 1.8) is ruled out with a high de¬ 
gree of confidence for our reference system, consistent with 
previous published results (B13,|Cappellari et al.|20T2 Con- 


roy van Dokkum||2Q12b Spiniello et al.||2014| . However, 

unlike previous results, we do not hnd any relation between 
the IMF slope and stellar velocity dispersion: dx/d\og{a) 
= for XLENS and dx/d\og{a) = -0.60tjj° for 

XLENS+S15. We note that the lack of a IMF-cr relation does 
not necessarily contradict previous studies, since this is the 
hrst time in which both the IMF slope and IMF low-mass 
cut-off are treated as free parameters and the range of stel¬ 
lar velocity dispersions is relatively small (280-380 kms“^). 
If we force a trend of x with cr, we must also have a trend of 
M\ow with cr. This indicates a degeneracy between the IMF 
slope and the IMF low mass cut-off, as expected. 

We hnd a marginally signihcant relation between IMF 
slope and effective radius: dx/d\og{Ref[) = O.GGIq'Jq for 
XLENS and dx/d\og{Res) = O.dTt^Jg for XLENS+S15. We 
conclude, albeit with low signihcance, that the effective ra¬ 
dius of ETGs might be a better indicator of IMF slope than 
their stellar velocity dispersion. Because in general the stel¬ 
lar dispersion correlates with effective radius, either can act 
as proxy for IMF slope, consistent with previous studies. 


Since both velocity dispersion and effective radius are 
global parameters, they might only be proxies for the physi¬ 
cal processes during star-formation. On the other hand, mass 
density might be a better indicator of the IMF slope and the 
low-mass cut-off. Figure shows the marginalized posteri¬ 
ors resulting from a lower-dimensional MCMC power-law 
model ht where we dehne the total dynamical density as 
p oc jR^^ in units of 1000 (kms“^)^/kpc^. 

We hnd that total density anti-correlates with the IMF 
slope, i.e. steeper slopes for less dense stellar systems, and 
that the lower mass limit increases as well for less-dense 
systems: dx/d\og{p) — —O.dTloJg for XLENS only and 


dx/d\og{p) = -0.19to\l for XLENS+S15 {dx/d\og{p) < 0 
at 87% CL). The relation is marginally signihcant, but when 
plotting the SSP-only inferred slopes of XLENS+S15 sys¬ 
tems against total density (Fig. 3), we see that the S15 lenses 
are all much denser and extend the anti-correlation trend al¬ 
ready seen from the XLENS systems. Moreover, dx/d\og{p) 
obtained with a simple least-squares ht of the data (solid 
lines, red for XLENS only and blue for the XLENS+S15 
case) agrees well with that obtained from the MCMC anal¬ 
ysis (dotted lines, same color code), but we emphasize that 
these hts are only a sanity-check of the MCMC results. 


4 SUMMARY DISCUSSION 

We combine the results of a state-of-the-art leasing and dy¬ 
namics analysis of nine massive lens galaxies of the XLENS 
Survey with corresponding inferences from a spectroscopic 
SSP study of line-strength indices focusing on the proper¬ 
ties of their IMFs. At the time of writing, this is the only 
method that allows the joint inference of the slope and low- 
mass cut-off of the IMF (see B13). 

We use a multi-variate, six-dimensional Markov Chain 
Monte Carlo analysis to derive the GL+D+SSP posterior 
PDFs on trends of the IMF slope and low-mass cut-off with 
stellar velocity dispersion, effective radius and total mass 
density. The main results of our analysis are as follows: 

(i) For a reference galaxy with a = 250 km s~^ and Reff = 

10 kpc, we infer X250 = 2.37to.i2 Mow, 250 = 0.13lto.o26 

M 0 , consistent with the results previously published in B13 
and with other published studies on massive ETGs. 

(ii) A marginally signihcant correlation is found with the 
size of the galaxies, i.e. systems with larger Reff prefer 
bottom-heavier IMFs. This correlation is similar to that with 
stellar velocity dispersion because more massive galaxies, 
statistically, also have larger effective radii. 

(iii) We hnd an anti-correlation between IMF slope and 
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Figure 3. IMF slopes of XLENS (black triangles) and S15 sys¬ 
tems (green squares) against total density. Solid lines show a least- 
square fit of the data, dotted lines show the relations obtained 
from the MCMC analysis. Red lines show the XLENS only case, 
blue lines show the XLENS+S15 case. The IMF slope plotted 
here is the one listed in Table. 1, obtained via SSP analysis alone. 


total (dynamical) density, defined as p = This trend 

could explain the seemingly contradictory result by |Smith^ 
Lucey k, Conroy| ( |2015 ) that suggests that some very high- 
velocity dispersion lenses have shallow IMFs slopes, counter 
to previous findings. Combining the XLENS+S15 sample, 
we find that the IMF slope anti-correlates (at 87% CL) with 
total mass density: dx/d\og{p) = —O.lQlQ ^g. 

The fully Bayesian analysis of the IMF parameters using the 
XLENS systems confirms our previous results in B13 that 
the IMF slope is inconsistent with that of the Milky Way at a 
very high confidence level for reference galaxies with a = 250 
kms“^ and 10 kpc. We find marginal (~ 1-cr) trends 

with effective radius and total mass-density. No trend with 
stellar velocity dispersion is found (although this is still con¬ 
sistent with our previous work). When we include the sample 
of S15, both trends become stronger; in particular the anti¬ 
correlation of the total mass density with IMF slope is found 
at the 87% CL. The S15 galaxies are much denser, with high 
stellar dispersion (cr = [356 zb 16,356 zb 18,320 zb 18]) but 
smaller effective radii — [6.6, 2.0, 6.0]), than the 

XLENS galaxies. If indeed this trend with mass density is 
genuine, it explains the seemingly discrepant results between 
most studies and that by SI5, but it also implies that a more 
physical (anti)correlation might exist between the low-mass 
IMF slope and the inner total mass density of ETGs. The 
reason for this anti-correlation remains an open question but 
it is likely related to the gas density during star formation. 


as suggested by recent theoretical papers (e.g., Krumholz 
2011 Hopkins||2012 Narayanan Da'^|2013 ). 
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